Doubly constrained bounds on the entanglement of formation 
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We derive bounds on the entanglement of formation of states of a 4 x N bipartite system using two 
entanglement monotones constructed from operational separability criteria. The bounds are used 
simultaneously as constraints on the entanglement of formation. One monotone is the negativity, 
which is based on the Peres positive-partial-transpose criterion. For the other, we formulate a 
monotone based on a separability criterion introduced by Breuer (H.-P. Breuer, e-print quant- 
ph/0605036). 
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PACS numbers: 03.67.Mn 

The nonclassical correlations of entangled quantum 
states have been of interest since the very inception 
of quantum mechanics 0, Q ■ Quantum information sci- 
ence has led to the idea that entanglement is a resource 
for information processing and other tasks. The abil- 
ity of quantum computers to solve classically hard prob- 
lems efficiently, the increased security of quantum cryp- 
tographic protocols, and the enhanced capacity of quan- 
tum channels all these arc attributed to entanglement. 
Investigating entanglement has led to new understand- 
ing of techniques such as the density-matrix renormal- 
ization group Q and of quantum phase transitions 0, 0] 
and properties of condensed systems jjj . Despite the im- 
portance of entanglement, however, characterizing and 
quantifying it in most physical systems that are of inter- 
est from an experimental standpoint remains a challenge. 

An important measure of entanglement for a pure 
state \^f) of two systems, A and B, is the entropy, 
—Tr (pa log pa), of the marginal density operator pa- We 
write this entropy sometimes as a function h(^>) and 
sometimes as the Shannon entropy H(fi) of the vector 
fi of Schmidt coefficients of This measure can be 

applied to bipartite mixed states by the convex-roof ex- 
tension of h(^>). The extended quantity, called the en- 
tanglement of formation (EOF), is defined as 



h(p) — min 

{pj.li"}} 



3 



(i) 

The EOF is a nonoperational measure of entanglement 
because the minimization over all pure-state decomposi- 
tions of p generally means there is no efficient procedure 
for calculating it. This minimization is the bottleneck in 
evaluating most nonoperational entanglement measures 
for mixed states. Consequently, bounding the EOF, in- 
stead of computing its value, becomes important. 

An alternate approach to quantifying entanglement is 
based on the use of positive (but not completely positive) 
maps on density operators [8j. A quantum state is sepa- 
rable if and only if it remains positive semidefinite under 
the action of any positive map. Given a positive map, 
we can construct a related entanglement monotone by 



considering the spectrum of density operators under the 
action of the map [jj 0] . Such monotones are typically 
much easier to calculate for general quantum states, be- 
cause they do not involve the convex-roof construction, 
and thus are said to be operational 

We can use the monotones constructed from positive 
maps and from other operational entanglement criteria as 
constraints to obtain bounds on nonoperational, convex- 
roof extended measures of entanglement. The complexity 
of the minimization in Eq. Q is reduced by solving it over 
a constrained set, instead of over all pure-state decom- 
positions. This was done in 0,0] for the EOF, using a 
single operational constraint. Our endeavor in this Let- 
ter is to carry this program forward. We first sketch a 
general scheme for many constraints, which we discuss 
further in , and then illustrate the general scheme for 
a particular case of two operational constraints. 

Let us say that /i , . . . , fx are operational entangle- 
ment monotones for a bipartite system. We gather 
their values for an arbitrary state p into a vector n = 
(m, . . . , Uk)- Their actions on pure states are func- 
tions of the Schmidt coefficients, i.e., fk(^) = Fk(t^) for 
k = 1,...,K. 

We are interested in a lower bound on the value of the 
EOF. Let us assume that for the state p, the optimal 
pure-state decomposition is p = Y^jPj\^){^\j gi ym g 
h(p) = J2j PjH(^) ■ Now define the function 

H(m)=min{lf(/x)|F fe (/x) = m fel k = 1, . . . , k} . (2) 

Notice that H(m) is defined only on the region of pos- 
sible values of m corresponding to pure states, a region 
we call the pure-state region. If H is not a monotoni- 
cally nondecreasing function of m, which we will call a 
monotonic function for brevity, we replace it with such a 
monotonic function if|(m), constructed by dividing the 
pure-state region into subregions on which subsets of the 
constraints are applied. We describe the procedure for 
constructing £zj(m) in detail in fl3| . 

Let H.(m) = co[i/j(m)] be the convex hull of 
iij(m), i.e., the largest convex function of K variables 



2 



(mi, . . . , m,K ) bounded from above by ffjfml We can 
show that H(m) is also a montonic function [13J, which 
can be extended naturally to a monotonic function on 
the entire space of values of m. Using Eq. J2J) and the 
convexity and monotonicity of Ti., we can write 



Kp) > 5>,W(n') > H[ ) > WW , (3) 



where we have used the convexity of the monotones fk 
to obtain J2jPj n k — nk - Knowing the easily calculated 
n for p thus leads to a bound on h(p). 

We now carry through the general program for 4 x N 
states using two operational entanglement monotones as 
constraints. Ours is the first instance of a doubly con- 
strained bound on an entanglement measure for a family 
of states. It g ives tighter bounds than those obtained 
previously [llj. 

The first monotone is the negativity [lOj , which is based 
on the Peres criterion 0] . The negativity of a bipartite 
state p is defined as Ut{p) = {\\p Ta \ \ — l)/2 where Ta is 
the partial transposition with respect to system A and 
the trace norm is defined as ||0|| = Tr(VOOt). For 
pure states, the negativity, in terms of the Schmidt coef- 
ficients, is given by nx = [Q2j \fWj) 2 ~ l]/2< 

We define a second monotone based on the <&-map in- 
troduced by Breuer 0. The action of the $-map on 
any state a is given by $(cr) = Tr(a)I — a — Va T V\ 
where the superscript T stands for transposition and V 
is a unitary matrix with matrix elements (j, m\V\j, m') = 
{— l) J ~ m <5 m ,-m' in the angular momentum basis {\j, m}}. 
The map $ provides, for any bipartite state p having a 
subsystem with even dimension greater than 4, a non- 
trivial condition for separability as (I ® < &)(p) > 0. 
The related entanglement monotone, which we call the 

-negativity, is defined for a general mixed state as 



n$(p) 



D(D-l) 



\\(i®mp)\\ 

D-2 



- 1 



(4) 



where D is the dimension of the smaller of the two sys- 
tems in the bipartite state p. The ^-negativity is a 
convex function of p. For 4 x N systems (N > 4), 
the ^-negativity for pure states, as a function of the 
four Schmidt coefficients, is n$ = 3^/ (pi + pa){p2 + Ma)- 
The ^-negativity for various states is given in |l3j . 

We can place bounds on the EOF of 4 x N states by 
using either n$ or nx as constraints. To find the bound 
with jit as the single constraint, which was done in [llj . 
one first finds the singly constrained function H(jit) of 
Eq. (J2J. This function being monotonic, but not convex, 
its convex roof gives the bound. For the 4 x N states we 
consider, the bound is given by 



H(n T ) 



H 2 ( 7 ) + (l-7)log 2 3, n T E [0,1], 
(nr- |)log 2 3 + 2 , n T € [1, §], 



(5) 
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FIG. 1: The shaded region is the pure-state region in the n$- 
riT plane for 4 x N pure states. The dashed lines are the 
monotone boundaries given by Eq. ® and by tit = n<s>/3. In 
the 2-constraint region between the monotone boundaries, we 
set -Hj(n) = H(n), and in the 1-constraint region above the 
upper monotone boundary, we set i?f(n) = H(ut)- 

where H2 is the binary entropy function and 7 = 
(V2n T + 1 + V9-6n T ) 2 /16. If instead we use n# as 
the single constraint, we first find the function H(n$), 
which being monotonic and convex, gives directly a dif- 
ferent bound on the EOF of 4 x TV states 0, 



Hfa) = H{n 9 ) = H 2 (a) 



1 + y/l - Anl/9 



1 and 3 v /(/ii -I- P4)(p 2 + M3) = 
P2, and P3 in terms of n$, nx, 



(6) 

We refer to 7i(n$) and 7i(nr) as singly constrained 
bounds on the EOF. We now proceed to place a doubly 
constrained bound on the EOF of 4 x N density operators 
by simultaneously using nx and n$ as constraints. 

Both n$ and nx take on values between and 3/2, so 
all 4 x N states lie in a square of side 3/2 in the n^-n^ 
plane. Not all points in the square correspond to pure 
states. Solving simultaneously the normalization con- 
straint Y)j—i Pj = 1 and the two constraint equations, 

Ej=iV^J = V2n T 4 
n$, lets us express p\ 
and p4. For some values of n$ and nr, there is no value 
for p4 for which the other three Schmidt coefficients are 
real numbers in the interval [0, 1]. 

To find the pure-state region, we look for the max- 
imum and minimum allowed values of nx for a fixed 
n$, assuming a pure state. To find the maximum, we 
apply the technique of Lagrange multipliers and obtain 
nx = 2rt$/3 + 1/2. The minimum lies on the boundary 
of allowed values of fx, with P3 = Pi = 0, and is given 
by nx = n$/3. The resulting pure-state region, shown in 
Fig. ^ is convex. The pure-state region is not convex in 
general, however; the subtleties this introduces into our 
program are addressed in . 

To find the doubly constrained bound on the EOF, 
we start with the function P|). specialized to our two 
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constraints, 

H(n<£, tit) = min i H(fi) 




V^7 = y/2n T + 1, 

3 

3 a/(a*i + Jh) (pi + ^3) = I • (7) 



It turns out that iJ(n$,riT) is not monotonic, so we 
must replace it with the monotonic function H^(n^,nT) 
discussed above. The procedure for constructing 
H-\(n$>,riT), depicted in Fig. ^ makes a connection to 
the singly constrained bounds. This connection is based 
on the fact that the minimum of any function subject to 
two constraints is greater than or equal to the minimum 
of the same function subject to only one of the^two con- 
straints. Thus we can say that H(n$,nT) > H(nx) for 
all n$ and H(n^,nr) > H(n<&) for all tit- 

The minimum of H(fjb) subject only to the tit con- 
straint, i.e., H{jit), occurs when the Schmidt coefficients 
are given by /x = (7, 7', 7', 7') ^lj with 7' = (1 — 7)/3. 
This corresponds to n$ = ^/2(2~f + 1)(1 — 7), thus defin- 
ing a curve in the n^-riT plane. Writing 7 in terms of Ut 
puts this curve in the form 



Along this curve, which we call a monotone bound- 
ary, the n$ constraint is automatically satisfied when 
H(fi) is minimized with respect just to the ut con- 
straint, which means that H(n^,nT) — H(nr) on this 
monotone boundary. To construct the required mono- 
tonic function, we set H^(n^,nx) — H(nr) when n$ < 
\/2(2"f + 1)(1 — 7), i.e., above this monotone boundary. 

Similarly the minimum of H(fi) subject just to the 
n$ constraint, i.e., H(n<s>), occurs when = (a, 1 — 
a, 0,0), which gives a lower monotone boundary ut = 
^ a(l — a) = n$/3. Along this line, the nx constraint is 
automatically satisfied when H([i) is minimized with re- 
spect just to the n$ constraint, which gives H (n$, n-r) = 
i?(n$) on this boundary. Since this lower monotone 
boundary coincides with the lower boundary of the pure- 
state region, it has no impact on defining i?|(n<j>, nx)- 

The definition of -ff-f(n) is depicted in Fig. ^ Be- 
tween the monotone boundaries, a region we call the 2- 
constraint region, we set (n) = H(n), and in the pure- 
state region above the upper monotone boundary which 
we call the 1-constraint region, we set -ffj(n) — H(tit)- 
The resulting function -ff-f(n) is monotonic throughout 
the pure-state region. 

We now focus on finding H(ja) in the 2-constraint re- 
gion. The method of Lagrange multipliers is not suit- 
able for finding the minimum Q because the problem is 
overconstrained. The equations obtained using Lagrange 



multipliers have a consistent solution only if n$ and ut 
are related as in Eq. JSJ, in which case H(n) = i?(ny). 
This does not mean that there is no minimum for H (fj.) 
for other values of n$ and ny, just that the minimum 
lies on a boundary of allowed values of ix. The bound- 
ary with three of the Schmidt coefficients being zero is 
the origin in the n$-% plane, where H(fi) = 0. The 
boundary with two zero Schmidt coefficients is the line 
Ut = n$/3, and along this line i?(n) = H(n$). 

The minimum of H(fi) in the remaining part of the 
2-constraint region can be found using a straightforward 
numerical procedure. As discussed above, the constraint 
equations can be solved to express \i\ , [ii, and /J3 in terms 
of n$, ut, and /J4. There are two distinct solutions, 
and . For a particular value of /J4, one or both of 
these solutions can be invalid in parts of the pure-state 
region because one or more of the three Schmidt coeffi- 
cients lies outside the interval [0,1]. For valid solutions 
we compute the entropy H(n). 
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FIG. 2: The part of the 2-constraint region covered by four 
values of /14. The two lines are the monotone boundaries. 

We first consider the boundary where one Schmidt co- 
efficient is zero by setting /X4 — in the solutions y^ x > 
and [j,( 2 \ Not all points in the 2-constraint region can be 
reached if we set 114 = 0. This is easily seen by noticing 
that the point n$ = — 3/2 corresponds uniquely to a 
maximally entangled 4 x N state, and for this state all 
four Schmidt coefficients have the value 1/4. Indeed, a 
continuum of points cannot be reached if we stay on the 
boundary defined by \i± — 0, so we increase the value of 
/i4 in small steps. The parts of the 2-constraint region 
that are covered by four values of /X4 are shown in Fig.0 

This numerical procedure gives us, for each point n = 
(n$,nT) in the pure-state region, the range of values of 
/J4 for which H(/j,^') and/or H(fi^) can be calculated 
at that point. The minimum of these entropies over the 
allowed range of values for \i± is the value of H(n). 

The function H(n) in the 2-constraint region is, as re- 
quired, a monotonic function of both n$ and tit- It is 
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FIG. 3: (Color online) The doubly constrained bound Ti.(n) 
on the EOF of 4 x N states. Also shown is a contour plot of 
the same function. 

extended to the monotonic function H-\ (n) on the the en- 
tire pure-state region using the procedure outlined above. 
The monotonic function Hj(n) is not convex, however, 
so we must compute its convex hull H(n). This can be 
done numerically, and it turns out that the difference be- 
tween H(n) and H-\(n) is quite small (~ 10~ 3 ), the two 
functions differing only in a small area near the max- 
imally entangled state. Had the pure-state region, on 
which -ffj(n) is defined, not been convex, 7i(n) would be 
defined on an extended convex domain 

To obtain a bound on the EOF of all 4 x N states, we 
have to extend H(n) outside the pure-state region to the 
rest of the n^-nx plane. The extension has to preserve 
the monotonicity of 7i(n) so that the string of inequali- 
ties in Eq. © holds. This is achieved by extending 7i(n) 
using surfaces that match the function at the lower and 
upper boundaries of the pure-state region. To preserve 
monotonicity, the surface added in the region below the 
lower boundary has zero slope along the nx direction, and 
the surface added in the region above the upper bound- 
ary has zero slope along the n$ direction. The resulting 
doubly constrained bound H(n) on the EOF is shown 
in Fig. [3J The figure indicates that the extension to 
the whole n$-ny plane produces a smooth and seamless 
surface. 

A third constraint based on the realignment crite- 
rion [lfil can be used to improve our bound on the 
EOF for certain classes of states. We can define the re- 



alignment negativity for a bipartite density operator p 
as n R = (\\K(p)\\ - l)/2, where [R-{p)] ijM = Pik,jl- For 
pure states, tir = hit- This means that in deriving the 
bounds, we could have redefined ut as max(rtx, tlr). 

In this Letter we focused on the derivation of a par- 
ticular doubly constrained bound on the EOF of 4 x 
N systems. Starting from the <J>-map introduced by 
Breuer [Til Il8l| , we defined an entanglement monotone, 
the <I>-negativity, and combined it with the usual negativ- 
ity to formulate a doubly constrained bound. We found 
that the pure-state region in the n^-riT plane is divided 
into sectors by monotone boundaries. The doubly con- 
strained pure-state marginal entropy is applicable only in 
the region between the monotone boundaries. In the re- 
maining portions of pure-state region, singly constrained 
entropies are applicable. Monotonicity and convexity dic- 
tate how to extend the bound to all states. We expect 
these features to persist for systems that are not 4xJV 
and for more than two constraints, in which case the 
monotone boundaries will generally be hypersurfaces. A 
sector in which an m-constrained marginal entropy holds 
will be bounded by sectors in which (m — l)-constraincd 
marginal entropies hold. These methods might provide a 
useful procedure for bounding the EOF and other convex- 
roof entanglement monotones. 
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